Semi— analytical solution of the Kondo model in a magnetic field 



C. Slezak, 1 - 2 S. Kern-em, 1 Th. Pruschke, 1 and M. Jarrell 2 

1 Center for electronic correlations and magnetism, Theoretical Physics III, 
Institute for Physics, University of Augsburg, 86135 Augsburg, Germany 
2 Department of Physics, University of Cincinnati, Cincinnati OH 45221, USA 

(Dated: February 1, 2008) 

The single impurity Kondo model at zero temperature in a magnetic field is solved by a semi- 
analytical approach based on the flow equation method. The resulting problem is shown to be 
equivalent to a resonant level model with a non-constant hybridization function. This nontrivial 
effective hybridization function encodes the quasiparticle interaction in the Kondo limit, while the 
magnetic field enters as the impurity orbital energy. The evaluation of static and dynamic quantities 
of the strong-coupling Kondo model becomes very simple in this effective model. We present results 
for thermodynamic quantities and the dynamical spin-structure factor and compare them with NRG 
calculations. 



I. INTRODUCTION 



The single impurity Kondo or s-d model (SIKM)0 
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is, together with its close relative, the single impurity 
Anderson model (SIAM)cl, one of the fundamental mod- 
els in the theory of correlated electron systems. It has 
been studied extensively over the past four decadesB, 
but despite its simplicity, no complete analytic solu- 
tion exists that provides information on both thermo- 
dynamic and dynamical quantities. For example, the 
Bethe ansatz solutionoa has solved all universal prop- 
erties of the Kondo problem including its high- to low- 
temperature crossover behavior. But the Bethe ansatz 
requires the limit of infinite conduction band width and 
cannot be used to calculate dynamical quantities beyond 
the low-energy limit. 

Thus, numerical methods, i4ik e Wilson's numerical 
renormalization group (NRG)M or quantum Monte- 
Carlo (QMC)u, in connection with maximum-entropy 
methodsQ have to be employed to access dynamical quan- 
tities on all energy scales. In both methods, evaluation of 
dynamical properties, and quantities related to them, like 
the Korringa-Shiba relation or the Friedel sum ruleo, suf- 
fer from unavoidable numerical errors. QMC simulations, 
in addition, cannot be used at zero temperature and are 
restricted to comparatively large values of JB Moreover, 
a reliable evaluation of single- and two-particle spectra 
and related quantities in an external magnetic field, as 
well as their comparison and interpretation within-aJocal 
Fermi liquid pictured become rather problematicDu, es- 
pecially in the limit of vanishing external magnetic field. 
Approximate analytical techniques, like perturbation ex- 
pansions or 1 /./V-expansionaj, typically only describe cer- 
tain properties of the Kondo model correctly. A notable 
exception in thisj-*espect is the so-called local moment 
approach (LMA)li! This perturbative approach is very 
accurate in most casesliilj-iBcluding those of nonvanishing 
external magnetic fields.BtL3 



Besides its relevance for the study of moment forma- 
tion in metals, the SIKM ([!]) has gained new importance 
as input for investigating non-dilute correlated electron 
systems like heavy fermion materials within dynamical 
mean- field theory (DMFT).El Here a reliable method 
for calculating single-particle correlation functions, espe- 
cially close to the Fermi energy, is extremely important. 

In this paper, we propose a new non-perturbative 
semi-analytical approach to the Kondo problem based 
on Wegner's flow equation methodllil and previous work 
on applications of flow equations to strong-coupling 
problemsE-3'EJl We will show that to a very good approxi- 
mation, many physical quantities of the Kondo model can 
be calculated from a resonant level model (RLM) , where 
the interacting features of the Kondo model are encoded 
in a non-constant effective hybridization function of this 
resonant level model. Surprisingly, this noninteracting 
effective model describes both universal low-energy prop- 
erties like the Wilson ratio as well as high-energy power 
laws and logarithmic corrections with very good accu- 
racy. Due to the noninteracting nature of this effective 
model this mapping allows immediate insights into the 
physics of the SIKM, for example the dependence of its 
static and dynamical quantities on a local magnetic field. 

After presenting our approach in the next section, we 
will discuss several static quantities at T = as a func- 
tion of a local magnetic field and derive analytical ex- 
pressions for their asymptotic behavior. As an example 
for a dynamical quantity, we will then discuss the spin- 
structure factor and the Korringa-Shiba relation. An 
outlook on potential future applications of our approach 
concludes this paper. 



II. MAPPING TO A RESONANT LEVEL 
MODEL 

A. Principles of the flow-equation method 

The general framework of the flow equation methodQ 
and its application to the Kondo model has been ex- 
plained in detail in Ref. flll Here we will only repeat the 
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main steps in order to make this paper self-contained, 
and refer to Ref. 16 for more details. 



The key idea of the flow equation approach consists in 
performing a continuous sequence of infinitesimal unitary 
transformations on a given Hamiltonian 



dH{B) 
dB 



= [r,(B),H(B)} 



(2) 



With an anti-Hermitian generator rj(B) the solution of 
equation (Q) describes a family of unitarily equivalent 
Hamiltonians H(B) parameterized by .the flow parame- 
ter B. By choosing r/(B) appropriatelytJ one can set up 
a framework that diagonalizes a many-particle Hamilto- 
nian H(B = 0), i.e. H(B = oo) becomes diagonal. 

The concrete realization of this approach for the Kondo 
model was discussed in Ref. The starting point is the 
bosonized formEZl of the Hamiltonian (|l|). Since we will 
be mainly interested in describing the basic ideas of our 
approach, we restrict ourselves to a linear dispersion re- 
lation. Notice, however, that the flow equation approach 
can also be used for a nontrivial conduction band den- 
sity of sta" 



as it does not rely on the integrability of 
the model.Ea With a linear dispersion relation the Kondo 
problem becomes effectively one-dimensional, the charge 
density excitations in (Q) decouple, and we only need to 
look at the spin density part 



H = Hn- 
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Here <j(p) = 



1 c„t — cIj_„i c„, ) are the bosonic spin den- 
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sity modes with the bosonic spin density field defined by 

®( x ) = ~i Eg^o — e~ lqx - a \ q \l 2 a(q). For simplicity we 
have set the Fermi velocity Dp = 1. a is proportional to 
the inverse conduction band width. All our latter results 
will be expressed as universal functions of the low-energy 
Kondo scale Tk, and we can consider (^) to be equivalent 
to our original Kondo Hamiltonian if Tk <C a -1 . 

Eq. (||) was used as the starting point H(B = 0) of 
the flow equation approach in Ref. |l6|. Away from the 
Toulouse point the unitary equivalence of the flow holds 
only approximately, but thiSpapproximation can be con- 
trolled by a small parameters and yields very accurate 
results. During the flow the Hamiltonian can be param- 
eterized as 

H(B) = H a +J2g P (B) (Cl(X(B))S- + h.c.) (4) 
p 

+ 5> p (5) [Cl(\(B p )),C p (\(B p ))} . 



Here B p = p~ 2 , and C^(A), C p (X) denote normalized 
vertex operators with scaling dimension A > in mo- 
mentum space, 
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that obey (0|C p (A) C\, (A)|f2) = 5 pp >6(p) and 

(0|Ct(A)C p ,(A)|n) = 5 pp ,9{-p). For the special 
case A = 1 they fulfill fermionic anticommutation 
relations {C|(l), C p ,(l)} = § pp i and can therefore be 
interpreted as creation and annihilation operators for 
fermions. 

In Ref .Jl6| the following flow equations for the param- 
eters in (|J)nave been derived 

dg P 
dB 
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and a differential equation for the flow of the scaling di- 
mension 



d\ 2 8ttA 2 (1 - A 2 ) v-^ , 
-dB= r(A 2 ) 1^9 q9 - q ka\' 



(7) 



It can be shownE3 that one always finds A — > 1 in 
the strong-coupling phase of the Kondo model, i.e. in 
the low-energy limit the vertex operators in (|J) become 
fermions. In the following we will use an improved version 
of the above flow equations by taking into account that 
all approximations should be performed with respect to 
the interacting ground state: It turns out that the only 
necessary modification in (||), @ and (Q) is that the ex- 
ponent in |ga| A " _1 gets replaced bx X 2 (B q ) — 1, i.e. it is 
not a running exponent anymoreO 



B. Equivalence to a resonant level model 

Now we will compare this system of differential equa- 
tions with the flow equations for a resonant level model 
(RLM). This will lead to the key result of this paper: the 
RLM can be used as an effective model for the compli- 
cated strong-coupling Kondo model. 

The Hamiltonian of the resonant level model is given 

by 
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Following the same flow equation approach as previously 
in the SIKM, we establish a solution to the RLM (|). A 
detailed description of the flow equation solution can be 
found in Ref. ^0|. One finds the following flow equations 
for the parameters in (pft 



dVk \ri ^2 

_ = -V k (e d -e k ) 



^ VuV 2 ^ + e q -2e d (g) 
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(10) 
(11) 



It should be noted that this yields the exact analytical 
solution. Having established the flow equations to solve 
both the SIKM and the RLM, respectively, one can now 
show an approximate equivalence of these two models. 
We introduce the substitution 
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and notice that with this substitution the two set of flow 
equations (|]j|) an d (PHll|) become equivalent for e d = 
in the RLM, with the exception of the logarithmic term in 
(|J) . Thus we now have established an approximate map- 
ping of the SIKM onto the RLM by means of (|l^) in the 
sense that their flow equation diagonalization is identical. 
We shall refer to this relation by introducing the effec- 
tive hybridization function A c ff(e) = 7rX)fc ^fc — e k) : 
the RLM with this non-trivial hybridization function 
can be used as an effective model for the SIKM in the 
Kondo limit (small coupling limit) PqJ — ► 0. Since this 
nonintcracting RLM is a simple, quadratic Hamiltonian, 
this mapping will allow us to read off and understand 
many properties of the complicated many-body Kondo 
physics in an intuitive and straightforward way. It will 
turn out that the deviations of A c ff(e) from a constant 
hybridization function encode the quasiparticle interac- 
tion and therefore the many-body Kondo physics in this 
quadratic effective Hamiltonian. 

Notice that the above mapping between the SIKM and 
the RLM becomes exact at the Toulouse pointed p a J = 
27r(2 — y/2) since X 2 (B) — 1 for all flow parameters B. 
One easily verifies that the effective RLM then has a 
constant hybridization function, A c ff(e) = jpgj 2 . In 
this case, our mapping just reduces to the observation 
already made by Toulouse that the partition function of 
the Kondo model for this specific coupling constant J is 
exactly equivalent to the partition function of a quadratic 
Hamiltonian Ell 

In order to specify the function A e ff(e) in the Kondo 
limit it is best to not directly use relation (|l^) , but to de- 
termine the effective hybridization function from match- 
ing a correlation function in the SIKM and the RLM. We 
have chosen the (S + (i)S _ (0))-correlation function, eval- 
uated it with respect to ([|) for B = 0, and then chose 
A e ff(e) in the RLM such that this coincided with the 
(d 1 '(t)(i(0))-correlation function. The resulting A c g(e) 
agrees with ( |l2"| ) in the high- and low-energy regimes, 
with deviations only in the crossover region. However, 
the mapping from the Kondo model to the effective RLM 
becomes better since this procedure manages to partly 
also take the logarithmic term in (^|) into account. 

The resulting effective hybridization function can be 
scaled into a dimensionless form with one dimensionful 



parameter A° ff oc Tk 

A eff (e) = A c ° ff A cff (e/A c ° ff ) 



(13) 



A c ff (x) is a universal function in the Kondo limit (paJ — > 
0). It is depicted in Fig. |l| for pgj — 0.1, and coincides 
with its universal form for |x| < 30 (|e| < 60Tk,ri|e. 
this should be sufficient for most practical purposesc3). 
For larger energies the effective hybridization function 
begins to cross over into linear behavior with logarithmic 
corrections depending on the bare coupling PqJ. The 
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FIG. 1: The dimensionless effective hybridization function 
A e ff(x) evaluated for p^J — 0.1 in the SIKM. The resulting 
function is symmetric and is only plotted for x > 0. It co- 
incides with the universal form for x < 30. The dashed line 
represents the fit ([L4|) and is nearly indistinguishable from the 
actual data. Notice especially the appearance of logarithmic 
behavior in the crossover region. 

following function provides an excellent fit (see Fig. |l|) 
A eff (.) = l + i ai ln(l+(£) 2 ) 



(14) 
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with the parameters from Table [l|. 



ao 


ai 


a 2 


0.829 


0.536 


0.00324 



TABLE I: Result of the fit (Q to the effective hybridization 
Ac«(x). 

A similar analysis based on the comparison of flow 
equations shows that the above mapping between the 
SIKM and the noninteracting RLM can be extended to 
the case of a Kondo Hamiltonian (Q) with a nonvanishing 
local magnetic field h 



-ffsiKM + gPB h S z 



(15) 
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by setting e d = gp B h in the RLM. However, the mapping 
with the above effective hybridization function becomes 
less accurate for <?//b|/j| ^ 7k due to the approximate na- 
ture of the flow equation solution We will discuss 
this point in more detail below. 

Summing up, as long as we are interested in static 
quantities in a local magnetic field smaller than approx- 
imately Tk and/or dynamical correlation functions for 
energies smaller than approximately 607k, we can use 
the RLM with the effective hybridization function (|l4| ) 
to describe the physics of the S1KM in the small coupling 
limit. The only undetermined parameter in the RLM is 
the energy scale Ag ff that explicitly depends on J. This 
overall energy scale is proportional to Tk- Notice that 
the non-perturbative behavior of this energy scale 

7k oc e" 1 /" ' 7 (16) 

follows correctly from the original flow equations (|| 0), 
compare Ref. |l6| 



that the low-energy excitations in the system are con- 
trolled by spin degrees of freedom, a well-known feature 
of the Kondo physics. However, in our approach this 
result can be read off directly from equ. (|17|). Second, 
Pd(0) oc 1/A e ff(0) oc 1/Tk, i.e. we obtain the correct 
scaling behavior for 7i mp directly from the behavior of 
A e fi(e). There is, however, a nontrivial correction com- 
ing from the factor in parenthesis in ([l7|). Note that for 
A e ff (e) =const. this correction is one, but for the strongly 
e-dependent A e ff (e) in Fig. |l| it is of the order of two. As 
we will demonstrate later, this difference is directly re- 
sponsible for obtaining the correct Wilson ratio in our 
approach. 

From the mapping S z = <vd — 1/2 it is easy to cal- 
culate Xzz(u + iS) = ~{gp B ) 2 {{d)d; <tfd)) ui+iS- Since the 
correlation function has to be evaluated within the RLM, 
one obtains for the imaginary part at T = 



x" z M = (3Mb) 2 t / duj' p d (J)p d (uj +J) 



(19) 



C. Calculation of physical quantities 

Once the mapping between the SIKM ([!]) and the effec- 
tive RLM (||) has been established, one can readily calcu- 
late physical quantities for the Kondo model. One com- 
plication arises from the fact that operators of the origi- 
nal SIKM have to be transformed by a unitary transfor- 
mation analogous to (||). In the language of the effective 
resonant level model they will thus in general correspond 
to more complicated many-particle operators. Since the 
intention of this paper is to demonstrate the potential of 
our mapping in a pedagogical setting, we will concentrate 
on two quantities that remain simple under these trans- 
formations: i) the z-component of the spin operator S z , 
which becomes S z = d) d — 1/2, and ii) the Hamiltonian 
itself. 

From the latter we obtain the internal energy f7i mp = 
(H — Ho) and the Sommerfeld coefficient, ji mp (h). A 
straightforward calculation in the noninteracting RLM 
yields 



7imp(ft) 



ir 2 B 



p d (0)(l-A'(0)) 



(17) 



where A'(w) denotes the derivative of 



A(w) = ip J de 



Aeff(e) 



and P J ... is the principal value integral. Here pd[t) is 
the impurity orbital density of states of the RLM 



A eff (e) 



7T {e-gush- A(e)) 2 + A c ye) 



(18) 



The result (|17j) for 

7imp has some interesting implica- 
tions. First, because <v is connected to 5+, it is apparent 



Again, this result provides direct access to an interpreta- 
tion of the behavior of \zz ( w + i&) in terms of the physics 
of the resonant level model. 



III. RESULTS 



One quantity that can be calculated analytically is the 
low-energy limit of the spin structure factor S(ui) == 



5(0) 



lim 



x2» 



(20) 



For a vanishing local magnetic field 5(0) is just the spin 
relaxation rate accessible in e.g. spin resonance experi- 
ments. With the result for Xzz( w ) from (fl9|) we obtain 



5(0) = (gp B ) 



1 



Ayo) 



tt ((gp B hr + Al a (0)Y 



(21) 



which leads to the curve shown in Fig. §. Eq. @ is of 
particular importance because it explicitly demonstrates 
universality, T k S(uj) — /(^/is^/Tk), and allows to di- 
rectly fit e.g. experimental data from ESR or NMR ex- 
periments and extract the Kondo temperature. Note fur- 
thermore that the result (|2l]) is not only valid in the 
Kondo limit, but also holds at the Toulouse point of 
the anisotropic Kondo model and everywhere in between. 
Since it does not depend on the details of A e g (e) it will 
also be true for general band structures eg in (|l]) and thus 
is eventually the result for 5(0) in DMFT calculations. 

The full frequency dependent x" z (w) has to be cal- 
culated numerically using the form of the effective hy- 
bridization function in Fig. [IJ. The results for three values 
of the external field, h = 0, gp B h = Tk and gp B h = 57k 
are displayed in Fig. ^|. These correlation functions pro- 
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FIG. 2: Universal curve for S(0) as function of the local mag 
netic field h. 
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FIG. 3: Xz Z (uj) for three characteristic local magnetic fields 
h — 0, gfJ.Bh — Tk and gptBh = 5Tk. 



vide a good example for the usefulness of our mapping 
to the effective RLM since one can directly interpret the 
structures and their frequency and field dependencies in 
terms of analytical formulas derived for the RLM. E.g. 
the high-frequency behavior of x" z ( w ) follows directly 
from equ. (|l9|) and the behavior of the effective hybridiza- 
tion function A e g(e) at large energies (which is linear 
with logarithmic corrections, see Fig. |l|): x"z( w ) decays 
like \/ui with logarithmic corpptions, in agreement with 
(expensive) numerical resultsci 

For the dependence of the dynamical susceptibility on 
the local magnetic field one makes use of the fact that 
the local magnetic field corresponds to the on-site en- 
ergy in the effective RML. Therefore it is obvious that 
the observed shift of the resonance peak in x"z( w ) IS due 
to the shifted center of the resonant level. Furthermore, 
the depletion of the maximum value is related to the 
decreasing occupation of the resonant level, which cor- 
responds directly to the increasing local magnetization 



in the SIKM. At the same time, one observes a decrease 
of the total spectral weight in x" 2 (w), which can be ac- 
counted for by a transfer to a finite expectation value of 
(S z ) in the SIKM. There is, however, also a non-trivial 
effect, namely the increasing broadening of the resonance 
peak with increasing magnetic field. For a RLM with a 
constant A G ff (e) such a behavior does not occur; it is en- 
tirely related to the fact that with increasing magnetic 
field the system starts to notice the energy dependence 
of the effective hybridization. 

The quantity not yet fixed in our calculation is Tk, 
or more precisely the proportionality constant in Tk cx 
A o ff(0). This can most conveniently be done by using 
Wilson's definition of the Kondo temperaturecl 



Xo (h = Q) = (g» B Y 



w 
47k" 



(22) 



where xo is the static magnetic susceptibility and w — 
0.413 the Wilson number, xo can be obtained from the 
imaginary part of the dynamic susceptibility (TL9) via 



Xo 



(23) 



and must in general be evaluated numerically. At the 
Toulouse point one can, however, give an analytic answer 
since A e ff (w) = const, and thus 



Xo(M = (5Mb ) 



1 



Aeff(O) 



(24) 



Therefcp 
relation^ 



IT W + A cff (0) 2 • 

at the Toulouse point the Korringa-Shiba 



R s = km 

27TXq LO 

is independent of the local magnetic field 



Rs 



(wb) 2 S(0) 



(25) 



(26) 



In the following we will discuss Xo(^) an d the 
Korringa-Shiba relation for the Kondo limit pqJ — > 0. 
The quantity Xo(^) is particularly convenient for a com- 
parison with NRG results. In Fig. El th e circles represent 
the values of xo(h) calculated via ( p3| ) with the effective 
hybridization function from Fig. |l|, and the full line rep- 
resents the result of an NRG calculation. We observe ex- 
cellent agreement for all values of the local magnetic field: 
notice that the curves agree without fit parameters. This 
example clearly demonstrates that the nontrivial form of 
the effective hybridization in Fig. [I] encodes the many- 
particle physics of the SIKM in a trivial noninteracting 
effective model. 

The result in Fig. ^ can readily be combinedj-svith re- 
lations ( p7| ) and (|2l]) to obtain the Wilson raticO 



Rw 



47T 2 fc| XimpW 
3(5Mb) 2 7imp(fr) 



(27) 
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FIG. 4: The magnetic susceptibility xo(h) from equ. ( p3| ) (cir- 
cles) and the same quantity obtained from an NRG calcula- 
tion. 



( |17| ) is very important to obtain this correct value for 
R-w(h = 0). Remarkably, our simple noninteracting ef- 
fective model therefore correctly describes the Wilson ra- 
tio in the Kondo limit (for not too large magnetic fields), 
which is a hallmark of strong-coupling Kondo physics. 

Let us finally analyze the accuracy of our effective 
model. Since Fig. ^ demonstrates that integral quantities 
like Xo(M are obtained with very good accuracy for all 
magnetic fields, one can infer from Fig. || that quantities 
depending on low-energy details in frequency space like 
7i mp and 5(0) are more susceptible to our approxima- 
tions for increasing magnetic fields. This suggests that 
for such low-energy quantities our effective model can be 
employed with very good accuracy (5% error) for mag- 
netic fields below 2k, and with good accuracy (20% error) 
still up to approximately 57k- 



IV. SUMMARY AND OUTLOOK 



Our results for the Wilson ratio and the Korringa-Shiba 
relation obtained within the effective RLM are collected 
in Fig. ||. For the Wilson ratio we woulcLactually have 
to calculate the quantity Ximp and not xo" However, for 
the case of small p$J considered here, both quantities 
are equivalent EB One observes that both i?w and Rs are 
independent of the local magnetic field up to approxi- 
mately g/J-sh 7k, and then start to decrease (Shiba 
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FIG. 5: Results for the Shiba ratio Rs (full line) and the Wil- 
son ratio Rw (dashed line) as a function of a local magnetic 
field. The correct limiting values at h — ► are missed by 
approximately 5%. 

ratio), respectively increase (Wilson ratio). The exact 
Bethe ansatz solutiona gives Rw(h) = 2 independent of 
the magnetic field strength (see also Ref. [12]), and lo- 
cal Fermi liquid theory yields Rs = 1 for h — > Ota Our 
limiting values as h — > miss these exact results by ap- 
proximately 5%. Notice that the term (1 — A'(0)) in 



Summing up, we have shown that the resonant level 
model with a nontrivial hybridization function A e g(e) 
can be used as an effective model for the single impurity 
Kondo model. The key observation was the fact that 
the flow equation solutions of both models are approx- 
imately identical if a suitable A c ff(e) is chosen for the 
RLM (Fig. 1). In this mapping the impurity orbital oc- 
cupation number — 1/2 plays the role of the Kondo im- 
purity spin 5 2 . The impurity orbital energy corresponds 
to the local magnetic field acting on S z . 

In contrast with the conventional approach where effec- 
tive models describe the vicinity of the low-energy renor- 
malization group fixed pointsuu, our effective model very 
accurately describes both certain low- and high-energy 
properties of the original Kondo model: compare for ex- 
ample our discussion of the dynamical spin-spin correla- 
tion function in Fig. |^. It also yields thermodynamic 
quantities that are in excellent agreement with much 
more expensive numerical methods (see Fig. |]). The 
nontrivial behavior of the effective hybridization function 
encodes the quasiparticle interaction, which leads to e.g. 
the correct Wilson ratio for small magnetic fields (with 
5% accuracy). Notice, however, that our effective model 
does not allow the correct evaluation of higher-order cor- 
relation functions beyond the low-energy limit, 

In conclusion, our approach describes many aspects 
of the complicated many-body Kondo physics for not 
too large magnetic fields within a simple noninteract- 
ing model. Therefore one can very easily and intuitively 
understand certain properties of the Kondo model, e.g. 
the dependence of correlation functions on a local mag- 
netic field (Fig. |^) . One main prospect of our approach is 
to look at other correlation functions using this effective 
model, in particular the T-matrix for applications in the 
framework of DMFT calculations. Future prospects also 
include cluster problems and the single impurity Ander- 
son model. Work along these lines is in progress. 
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